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An algebraic analysis of Grover's quantum search algorithm is presented for the case in which 
the initial state is an arbitrary pure quantum state \tp) of n qubits. This approach reveals the 
geometrical structure of the quantum search process, which turns out to be confined to a four- 
dimensional subspace of the Hilbert space. It unifies and generalizes earlier results on the time 
evolution of the amplitudes during the search, the optimal number of iterations and the success 
probability. Furthermore, it enables a direct generalization to the case in which the initial state is 
a mixed state, providing an exact formula for the success probability. 
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Oh Grover's quantum search algorithm 0, Q exemplifies the potential speed-up offered by quantum computers. It 
| also provides a laboratory for the analysis of quantum algorithms and their implementation. The problem addressed 
by Grover's algorithm can be viewed as trying to find a marked element in an unsorted database of size N. While 
a classical computer would need, on average, N/2 database queries (and N queries in the worst case) to solve this 
problem, a quantum computer using Grover's algorithm, would accomplish the same task using merely 0(^/N) queries. 
This proves the enhanced power of quantum computers compared to classical ones for a whole class of oracle-based 
problems, for which the bound on the efficiency of classical algorithms is known. Moreover, it was shown that 
Grover's algorithm is as efficient as theoretically possible I4L A variety of applications were developed, in which the 
algorithm is used in the solution of other problems @, 0, UH IH EH Ej) • 

Several generalizations of Grover's orig inal algorithm have been developed. The case in which there are several 
marked states was studied in Rcfs. was shown that when there are r marked states, Grover's algorithm can 

lO ' find one of them after T = O(^jNjr) queries. A further generalization was obtained by allowing the replacement of the 
Hadamard transform, used in the original setting, by an arbitrary (but constant) unitary transformati on |l4lll5l lT^ | . 
as well as by the replacement of the it inversion by an arbitrary (but constant) phase rotation ^3 Il8| ^ Another 
O ,' generalization was obtained by allowing the replacement of the uniform superpo sition of all basis states, used as the 
initial state of the algorithm in the original setting, by an arbitrary pure 19, 20] or mixed [21| quantum state. It was 
shown that the optimal time to perform the measurement that concludes the operation of the algorithm is independent 
of the initial state. However, the probability of success, P s , is reduced, and its value depends on the initial state. An 
^ | explicit expression for P s in terms of the amplitudes of the initial state was found |22| . This generalization provides 
an operational measure of entanglement of pure multi-partite quantum states [2^. 12 ll l25l] . 

In this paper we introduce an algebraic approach to the analysis of Grover's quantum search algorithm with an 
arbitrary initial quantum state. This approach reveals the geometrical structure of the search process, which turns 
out to be confined to a four dimensional subspace of the Hilbert space. This approach unifies and generalizes earlier 
results on the time evolution of the amplitudes during the search, the optimal number of iterations and the success 
probability. Furthermore, it enables a direct generalization to the case in which the initial state is a mixed state, 
providing an exact formula for the success probability. 

The paper is organized as follows. In Sec. II we briefly describe the algorithm. The algebraic analysis is presented 
in Sec. Ill for the general case that involves several marked states with an arbitrary pure state as the initial state. 
Special cases such as the case of a single marked state are considered in Sec. IV. The generalization to mixed initial 
states is presented in Sec. V. The results are summarized in Sec. VI. The detailed calculation of the success probability 
of the algorithm is given in the Appendix. 



II. THE QUANTUM SEARCH ALGORITHM 

Consider a search space D containing N elements. We assume, for convenience, that N = 2™, where n is an integer. 
The elements of D are represented using an n-qubit register containing the indices, i — 0, ...,N — 1. We assume 
that a subset of r elements in the search space are marked, that is, they are solutions of the search problem. The 
distinction between the marked and unmarked elements can be expressed by a suitable function, / : D — > {0, 1}, 
such that / — 1 for the marked elements, and / — for the rest. The search for a marked element now becomes a 
search for an element for which / = 1. To solve this problem on a classical computer one needs to evaluate / for each 
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element, one by one, until a marked state is found. Thus, on average, N/2 evaluations of / are required and N in the 
worst case. For a quantum computer, on which / to be evaluated coherently, it was shown that a sequence of unitary 
operations called Grover's algorithm can locate a marked element using only 0(^/N/r) coherent queries of / 0,0- 

To describe the operation of the quantum search algorithm we first introduce a register, \i) = \i\ . . . i n ), of n qubits, 
and an ancilla qubit, \q), to be used in the computation. We also introduce a quantum oracle, a unitary operator 
O, which has the ability to recognize solutions to the search problem. The oracle performs the following unitary 
operation on computational basis states of the register, \i), and the ancilla, \q): 



0\i)\q) = \i)\q®f(i}) (1) 

where ® denotes addition modulo 2. The oracle recognizes marked states in the sense that if |i) is a marked element of 
the search space, namely f(i) = 1, the oracle flips the ancilla qubit from |0) to |1) and vice versa, while for unmarked 
states the ancilla is unchanged. The ancilla qubit is initially set to the state \—) q = (|0) — |l))/\/2. With this choice, 
the action of the oracle is 0\i)\—) q = (— ) q . Thus, the only effect of the oracle is to apply a phase of —1 if 

x is a marked basis state, and no phase change if x is unmarked. The state of the ancilla does not change. 

Grover's search algorithm may be described as follows: Given an oracle O, whose action is defined by Eq. and 
n + 1 qubits in the state |0)®"|0) g , the following procedure is performed: 

1. Initialization: Apply a Hadamard gate W = A= ( } } 1 ) to each qubit in the register, and the gate WX to 

the ancilla, where X — ( ? J ) is the not gate, and we write matrices with respect to the computational basis 
(|0), |1)). The resulting state is \rj)\— ) q , where 

1 N ^ 

I,) = ^ £ ». (2) 



2. Grover Iterations: Repeat the following operation r times (where r is given below). 

(a) Apply the oracle, which has the effect of rotating the marked states by a phase of tt radians. Since the 
ancilla is always in the state \—} q the effect of this operation may be described by the unitary operator 
Im = I — ^XimeM \m)(m\, acting only on the register, where I is the identity operator. 

(b) (i) apply the Hadamard gate on each qubit in the register; (ii) apply the operator Iq = I — 2 1 0) (0 1 which 
rotates the 1 00 . . . 0) state of the register by a phase of tt radians, (iii) Apply the Hadamard gate again on 
each qubit in the register. 

The resulting operation is —WIqW = — I + 2 1 77) (77 1. When this operator is applied on the state J^i ^ 
results in the state — ai)|z), where a = ai/N. Thus, each amplitude is rotated by tt around the 

average of all amplitides of the quantum state. 

3. Measure the register in the computational basis. 

The combined operation on the register in one Grover iteration is given by Q — — WIqWIm- The optimal number 
of iterations before the measurement is 




where [a^J is the largest integer which is smaller than x [2I I3L Il3j. Moreover, at this optimal time a marked state 

can be found with almost certainty, or more precisely with probability P s = 1 — O (l/VN\ . With this performance, 

Grover's algorithm was found to be optimal in the sense that it is as efficient as theoretically possible Q. Note that the 
probability P s ~ 1 can be achieved only for specific initial states such as the one produced in step 1 of the algorithm 
above. If this initial state is replaced by an arbitrary quantum state, the probability of succes, P s , is reduced [lflll^ . 

The time evolution of the amplitudes of the marked and unmarked states during Grover's iterations was studied 
in Ref. |20| for an arbitrary pure initial state In particular, the optimal number of iterations and the success 
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probability were calculated, and found to depend on the specific choice of the set of marked states. Obviously, in a 
search process, the marked states are not known. Thus, the success probability should be averaged over all possible 
choices of the set of r marked states [22j . The results are that for any initial state the optimal number of iterations 
is given by Eq. J3J). The success probability, P s can be expressed in terms of the average of all amplitudes of the state 
\ip). In the next Section we introduce the algebraic approach to the analysis of Grover's search algorithm. 



III. ALGEBRAIC ANALYSIS OF THE QUANTUM SEARCH PROCESS 



Consider a search using Grover's algorithm, for one of r marked states in a space of N = 2" computational basis 
states, where n is the number of qubits in the register. The initial state is 



N-l 

= E a M ( 4 ) 

where is the amplitude of the basis state \i). Denote the set of indices of the marked states by Ai. The amplitudes 
of the marked states will thus be a mi m £ M. The complementary set, of unmarked states is denoted by Li, namely 
the amplitudes of the unmarked states are a u , u S Li. Thus, for a given choice of the set of marked states, the initial 
state can be expressed by 



E a m \m) + 2J a u \u). (5) 



The amplitudes a m and a u satisfy: 



Ei 



= 1. 



(6) 



Their averages are given by 



a>M = - E a ™ ( 7 ) 



for the marked states, and 



a u = 77 E Qu ( 8 ) 

N — r 



for the unmarked states. 



A. Construction of a four-dimensional subspace 

The initial state \tp) can be projected onto the subspaces spanned by the marked and unmarked basis states, giving 
rise to the normalized projections 



\<P 



Ml 



IP* 



E< 



(9) 



and 



(10) 
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respectively, where 

Po= i a ™i 2 - ( n ) 

The state can now be written in the form 



IV) = VPolM + Vi-Wtf), (12) 

where |0m) and \<pjj) are orthogonal to each other. Similarly, the equal superposition state can be expressed in 
the form 



where 



\v) = \I^\Vm) + J1-^\vu), (13) 



= 7f E M ( 14 ) 



and 



From now on we will refer to the plane spanned by \t)m) and as the Grover plane. Using the Gram-Schmidt 

procedure we extract from \4>m) and \<j>u) two new states \tpAt) and \4>u), that are perpendicular to \t]m) and |»7t/). 
These states are given by 

\Vm) = ; (16) 



\/i - I(^m|<Am)| 5 



and 



= ^ PES (i?) 



Using the fact that 



T _ 



Pi 







N - r _ 



we can write \4>m) and \ipu) more explicitly as 

(-/^I^m) - Vra M \r]M)) 



\i>M) 



\JPq- r\a M \' 



\*l>u) = - 1 — -■ (19) 

'l-P -(N -r)\au\ 2 
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Thus, Eq. I)12|l now takes the form 



\J 1 Pq- r\a M \ 2 \ipM) + \fl - Po - (N - r)\au\ 2 \ipu) + VW^rau^m) + v^mI^m), 



(20) 



where the state vectors \?Pm), \ipu), \vu) and \t)m) define a set of four orthonormal basis vectors. In particular, \ipM) 
and \t]m} span the subspace of marked states, while \ipu) and \i]u) span the subspace of unmarked states. Also, IV'm) 
and \ipu) a re perpendicular to I77). We will show below that for any initial state iterations of the quantum search 
always preserve the subspace spanned by these four vectors. 



B. The time evolution of the state vector 

Denoting the Grover iteration by Q, the state of the register after t iterations is 

\g(t)) = Q t m- 

Using the structure of the operator Q presented above, we obtain that for any state 



■ 2 ( <t#> - 2 £ (r,\m) (m|V) j \v) + 2 ]T { 



Applying Q on the vectors that span the four-dimensional subspace we obtain 



Q\vu) 
Q\vm) 



\iPm) 



l-^)l'/.u). 



(21) 



(22) 



(23) 



The Grover iteration Q acts as a linear transformation within a four dimensional vector space, spanned by \iI>m)i 
\ipu), \Vu) and \t]m}- For the analysis presented below it is convenient to use the vector representation 



IV 



lit 



\Vm) 



(24) 



in which, according to Eq. lj2U|) 



The matrix representation of Q is: 



( 



- r\a M \ 
I - Po - {N - r)\auf 

yN — rau 
Jra M 



(25) 



Q = 



100 

0-10 

cos uj — sin oj 

smut cos uj 



(26) 
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where 



and 



2r 

JV 



si,)u -' = 2t/ ^( 1- ^ 



(27) 



(28) 



Therefore, 



1 



Q 



Q 





(-1)* 
cos(wt) 
y sin(cjt) 






— sin(cjt) 
cos(wi) y 



(29) 



and the state of the register after t iterations is 



!<?(*)) 



■\JPo-r\a M \ 2 \i> M ) + (-l)Vl-fl)-(JV-r)|au| a |Vtr) 
+ K/ JV — raj/ cos(wi) — \frau sin(wi)^ j^;/) 

+ JV" — raj/ sin(wt) + ■sfrau cos(wi)^ |jjm)- 



(30) 



The four dimensional space in which the dynamics is confined includes the 'plane of marked states, spanned by IV'm) 
and \t]m) as well as the plane of unmarked states spanned by \tpu) and \f}u)- 

C. The success probability 

The probability of success P s (t) of a measurement taken after t iterations is given by the projection of \g{t)} on the 
plane of marked states: 



P s (t) = ^0 + 7. 



(JV~r)|%| 2 -r|a M | 2 



{N -r)\a v \ 2 -r\a M \ 2 



cos(2ujt) 



+ \ V r ( N - r ) (o-uO-m + a* M au) sin(2wt). 



(31) 



From Eq. I|3U|) it is easy to see that \g(t)) exhibits a rotation with angular velocity u> in the Grover plane spanned by 
\t}u) and \t]m}- This rotation causes an exchange of probability between the planes of marked and unmarked states. 
Thus, the success probability is limited by the projection of the initial state \if>) on the Grover plane. The projection 
of \ip) along \tpu)i which is perpendicular to the Grover plane, represents a lost probability that cannot be transformed 
into the plane of the marked states. This provides an upper bound on the success probability P s (t), of the form 



P ma x = Po + (N-r)\au\' 
Similarly, a lower bound on the success probability 



(32) 



J^min = Pa - r\a M \ 



(33) 
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can be obtained using the projection in the direction of \iPm)- Furthermore, by using the explicit forms of the states 
that span the four dimensional subspace and Eq. I|3U[) . the following expression for the temporal evolution of the state 
vector is obtained 

!<?(*)) = E (*(*) + Aa "0 l m ) + E ( f W + (-!) tA ««) l«>- (34) 



m£LA4 u£U 



The variables 



N-r _ 

~ \ a u sin(o;t) + a« cos(wi) (35) 



and 



lit) = ay cosfwt) — 4 / clm sm(ojt) (36) 

f iV-r 

are the average amplitudes, after t iterations, of the marked and unmarked states, respectively. Also, 

Aa m = a m - clm, m e M 

Aa u — a u — au, u<ElA (37) 

represent the initial deviations of the amplitudes a m and a u from the averages of the marked and unmarked amplitudes, 
respectively poj . 

The analysis presented above applies to any specific choice of the set of marked state. In practice, the set of marked 
states is not known. Thus, the success probability of an actual search is the average of P s {t) over all possible choices 
of the set of r marked states. In the limit of a large search space, where r N, it takes the form (see Appendix): 

(P s (t)) = N\a\ 2 sin 2 [uj(t + 1/2)] + O (r/N) , (38) 

where 

a = — E a< ( 39 ) 

i=0 

is the average over all the amplitudes of the initial state \tp). The optimal number of iterations, r, is given by Eq. J3J|, 
namely, is identical to that obtained for the original algorithm. The probability of success of a measurement taken 
after r iterations is 

P max = N\a\ 2 + 0(r/N). (40) 
It can also be expressed by -P max — \(-q\ip)\ 2 + 0{r / N) p^ . 

IV. ANALYSIS OF SPECIAL CASES 



Below we consider special cases in which the analysis can be simplified and the dimension of the subspace in which 
the Grover iterations operate is reduced. 
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A. Quantum search with a single marked state 

Consider the quantum search algorithm with an initial state and a single marked state |m). In this case the 
subspace of the marked states is one dimensional. Therefore, the normalized projection of the equal superposition 
state | rf) on the subspace of marked states is given by 

\vm) = \m). (41) 

Since the one dimensional subspace of the marked states does not include any perpendicular direction \iPm) = 0. The 
average amplitude of the marked states is clm = »m, where a m is the amplitude of the marked state \m) in the initial 
state \ip). Thus, the state \g(t)) of the register after t Grover iterations [Eq. (|30|l ] takes the form 

\g(t)) = (VN — lay cos(ujt) — a m sm(ujt)j \riu) + (y/N — lau sin(cji) + a m cos(wt)^ |m) 

+ (-l)V 1 -( iV - 1 )l^| 2 -| a ™| 2 \^u)- (42) 

Clearly, \g(t)) is confined to the three dimensional subspace spanned by \f]u)i \m) & n d \ipu)- 

In case that all the amplitudes are real, the dynamics of the quantum search can be viewed as a rotation of the 
state vector \g(t)} around the bases of a cylinder, jumping from one base to the other at each time step. The axis 
of the cylinder is in the direction of \ipu)- Its bases are in the Grover plane, namely spanned by \rju) and |m). The 
radius of the cylinder is given by the projection of \g(t)) on the Grover plane: 



R = \J(N -l)\au\ 2 + \a m \ 2 . (43) 
The length of the cylinder is twice the size of the projection of | <?(£)) in the direction of IVv), namely 



L = 2y/l - R 2 . (44) 

The dynamics of the quantum search consists of rotations of the state vector \g(t)) around the \ipu) axis, at an angular 
velocity lo, combined with switching positions between the two bases. At even time steps \g{t)) points towards the 
upper base, while at odd times it points towards the lower base. The optimal measurement time is when the vector 
\g(t)) is exactly above (or below) the \m) axis. 



B. Quantum search with an initial state in the Grover plane 



Consider a quantum search with a certain set of r marked states, such that the initial state \ip) is in the Grover 
plane. In this case \ip) can be represented by 



\1>) = a\riu) + P\Vm) (45) 

where a and (3 are complex amplitudes that satisfy \a\ 2 + |/3| 2 = 1. Note that \ip) has no component perpendicular 
to the Grover plane either in the subspace of marked states or in the subspace of unmarked states. Therefore, 
IV'm) = iV't/) = 0. The average amplitudes of the marked and unmarked basis states in are au = P/y/r and 
ajj = a/y/N — r, respectively. The state vector after t iterations takes the form 



\g(t)) — [acos(uit) — /3sin(cj£)] \i]u) + {asin(ut) + /3cos(ujt)] \t]m), (46) 

namely, it is confined to the Grover plane, where it rotates with angular velocity lo. The success probability P s (t) 
oscillates periodically between P m ; n = and P max = T namely the algorithm is optimal. For example, the equal 
superposition state \rf) is a special case within this category, where a = \J N — r/N and = y/r/N . 
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C. Quantum search with an initial state perpendicular to the Grover plane 

Consider an initial state which is perpendicular to the Grover plane, namely a state that satisfies (t)m\'>P) = 
(Vu\ip) — 0- In this case Hm = ajj = and 

|fl(t)) = a/Wa/) + (-!)* VT^PoHu). (47) 

This state vector is confined to a cycle of period 2, namely \g(t + 2)) = | <?(£)). Recall that Grover's algorithm generates 
a flow of probability from the subspace of unmarked states toward the subspace of marked states only within the 
Grover plane. In this case \ip) is perpendicular to the Grover plane. Therefore, the probability of success remains 
unchanged, namely P s (t) = Po, making the quantum search process useless. Two special cases can be identified. In 
case that du = and a m = for all m £ A/, the success probability P s {t) = at any time t. In the opposite case 
where am = and a u = for all u £ U, P s (t) = 1 at all times. 

V. QUANTUM SEARCH USING MIXED STATES 

The quantum search with an initial state which is a mixed state was studied before for a specific choice of the set 
of marked states j2]| . Here we extend the analysis to the general case in which the set of marked states is unknown 
and calculate the success probability. We first analyze the case of a genral mixed states and then discuss some special 
cases. 



A. General analysis for arbitrary mixed states 

Consider a quantum search using the mixed state 

A) = X>/^>^Vl (48) 

of n qubits as the initial state of the register. In this representation is the probability that corresponds to the pure 
state j^) in the ensemble and X^fV = ^ ne P ure states take the form 

|VV> = E <VI J >- ( 49 ) 

i=0 

The average amplitude in the state is 



The density operator after t iterations is given by 

pit) = = J2pMt)){gM (51) 

where |<7^(i)) = Q \ip^. The measurement of the register in the computational basis is represented by the operator: 



N-l 

M = Mi, 

i=Q 



(52) 
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where Mj = is the measurement operator that corresponds to the outcome i. The success probability of a 

measurement taken after t iterations is 

P s (t) = ]T Tr (MlM m p{t)) = 5^p M P.(^,t) (53) 

mEM H 

where 

P,fcM) = J2 IH5m(*))| 2 . (54) 
Averaging the probability of success over all possible choices of the set of r marked states, where r <C N, gives 

(P s (t)) = 7V|a| 2 sin 2 [w(t + 1/2)] + 0(r/N), (55) 

where 

|a| 2 = ^p M |a M | 2 . (56) 

Thus, the maximal success probability 

Pmax = N\a\ 2 + 0(r/N) (57) 

is simply the weighted average of the success probabilities obtained using the pure states as initial states. Using 

Eq. (|56J) and the fact that | (77 1 V'm) 1 2 = ^I^m| 2 we nnc ^ that Pnax — ('ylPol 7 ?)- This is, in fact, the square of the fidelity 
of po and \rj), namely, 

Pn^ = F 2 (\ v ),p ). (58) 

B. Search using only part of the qubits in the register 

Consider a quantum register of n + k qubits, which is in a pure state \ip). Suppose that the qubits are divided 
between two parties such that Alice gets the hrst n qubits and Bob gets the rest of them. The state of the register 
can be expressed by 

K-l N-l 

W=EE b^A^B (59) 

fj,=0 i=0 

where the subsystem of Alice is spanned by the basis \i)a, i = 0, 1, . . . , N — 1, where N = 2™, while the subsystem of 
Bob is spanned by \p)b, M = 0,1,.. ■ ,K — 1, where K = 2 . The normalization condition is 

K-l N-l 

EEi^i 2 = i (60) 

fj,=0 i=0 

and 

K-iN-l 

^M^^^ (61) 

M =0 i=0 
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is the average amplitude. Eq. I|59(l can be written in the form 



where 



K-l 

c Ma\»)b (62) 

fj,=0 



K-l 

E M 2 = L (63) 
The pure state \ipfi) a of party A, that corresponds to the computational basis state \[i)b of party B is given by 

N-i 

\^^)A=Y. a ^\j)A (64) 
j=o 

where b^t = c^a^. This state satisfies the normalization condition 

JV-l 

|2 



£M 2 = 1 (65) 



i=0 

and its average amplitude is 



1 



TV 

i=0 



The measurement statistics for Alice is given by the reduced density operator: 

K-l 

p A = Tr B |V)<^| = E pJ^A a(^\ (67) 

where Tr^ is a trace over the state of Bob and = |c M | 2 . 

Consider a quantum search in a search space of size NK using the entire register of n + k qubits in the state \ip). 
The success probability is 



= NK\W = ^ 



K-l 1 
fj,=0 



O 



Now, consider the case in which Alice performs a quantum search in a space of size N using her n qubits in the state 
p A . The success probability of Alice's search is 

K-l K-l 

PL* = «E PM 2 = N E \c^\ 2 + 0(r/N). (69) 

Using the inequality |ir| 2 < \x\ 2 , for a random variable Xk, where equality is obtained only if all the Xk's are equal, we 
find that > -P^ax- This means that for a register in an entangled state reducing the search space increases 

the success probability. Equality is obtained only if there is no entanglement between the systems held by Alice and 
Bob, and Bob's system is in the equal superposition state \t))b- Thus, one can always add or remove unentangled 
qubits in the state (|0) + |l))/v2 without changing P ma x- However, adding entangled qubits reduces P max while 
removing such qubits increases it. 
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C. Search using an initial pseudo-pure state 

Consider a search using an initial state which is a pseudo-pure state of n qubits 

p e = (l-e)^ + eM<V|. (70) 

where / is the TV-dimensional unit matrix representing the maximally mixed state, is a pure state and < e < 1. 
The maximally mixed state can be represented by 

N-l 
i=0 

while the pure state = a,i\i), has an average amplitude a^,. Using Eq. I|55|) we find that 

P max = ^+^VM 2 + 0{er/N). (72) 

Therefore, the success probability is simply reduced by a factor of e vs. its value in the case that the initial state is 
the pure state \4>). 

VI. SUMMARY 

We have introduced an algebraic approach to the analysis of Grover's quantum search algorithm with an arbitrary 
initial quantum state. This approach reveals the geometrical structure of the search space, which turns out to be a four 
dimensional subspace of the Hilbert space. This approach unifies and generalizes earlier results on the time evolution 
of the amplitudes during the search, the optimal number of iterations and the success probability. Furthermore, it 
enables a direct generalization to the case in which the initial state is a mixed state, providing an exact formula for 
the success probability. 



APPENDIX A: THE AVERAGE SUCCESS PROBABILITY 



The probability of success of the quantum search algorithm after t iterations, for a given choice of the set of r marked 
states is given by Eq. Since the set of marked states is not known (although their number, r, is specified) the 

actual success probability is the average (P s (t)) over all possible choices of the set of marked states: 



(Ps(t)) 



1 r 

2 



(N-r)(\au\ d 



(N-r)(\au\ 



r{\aM\ 



r{\a,M 

2 



(Po 

i r 

~2 

+ 2 V r ( N - r ) ((o-uO'M) + {a*uau)) sin(2wi). 



cos(2wi) 



(Al) 



We will now evaluate the averages (Po), (|a[/| 2 ), ( | a a# | 2 ) and (o^cim)- The average a over all amplitudes in \ip), 
satisfies the inequality \a\ < 1/N where equality is obtained only for the equal superposition state |?7). The average 
of Po over all possible choices of the r marked states is 



^AT-l N-l 



mEM r j=0 

where the binomial coefficient 



(Po) = < E M 2 > = -gr E M 2 = jj> ( A2 ) 
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is the number of different sets of r objects that can be picked out of N distinguishable objects. The second moments 
of the amplitude distribution are 

(l%| 2 > = 7—^2 ( E a «i E a U2 *) = \a\ 2 + 0(l/N?) 
(\a M \ 2 ) = \a\ 2 (l-l)+^ + 0(l/N*) 

(%«M> 

v ; meM ueu 

(a M a*u) = (a^a A f)* = |a| 2 + 0(l/7V 2 ). (A4) 



^rj^-A E a ™ E o - n 2 + (w 2 ) 



Substitution of the above averages in Eq. (|A1|) yields 

(P s (t)) = 7V|a| 2 sin 2 [c;(i + 1/2)] + ^ (l - 7V|a| 2 ) + (l/N) , (A5) 
In the limit of a large search space and r <C N, the expression for the success probability is reduced to 

(P s (t)) = N\a\ 2 sm 2 [Lu{t + 1/2)] + O (r/N) . (A6) 
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